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Abstract 



We prove a global in time existence and uniqueness theorem for the 
in : initial value problem for the relativistic Boltzmann equation on the flat 

^Tj ■ Robertson- Walker space-time, in a functional framework appropriate to 

the coupling with Einstein's equations. We had nowhere to restrict the 
size of the initial data which can hence be taken arbitrarily large. 



1 Introduction 



The Boltzmann equation is one of the basic equations of the kinetic theory. 
This equation rules the dynamics of one kind of particles subject to mutual 
collisions, by determining their distribution function, which is a non-negative 
real-valued function of both the position and the momentum of the particles. 
The distribution function is physically interpreted as the "probability of the 
jj] ■ presence density" of the particles in a giving volume during their collisional 

evolution. 

In the case of binary and elastic collisions, the distribution function is de- 
termined by the Boltzmann equation through a non-linear operator called the 
"Collisions operator" , that describes, at each point where two particles collide 
with each other, the effects of the behaviour imposed by the collision to the 
distribution function, taking into account the fact that the momentum of each 
particle is not the same, before and after the collision, only the sum of their two 
momenta being preserved. 

Due to its importance in the kinetic theory, several authors studied and 
proved local and global existence theorems for the Boltzmann equation, in both 
the non-relativistic case, that considers particles with low velocities, and in the 
relativistic case, that includes the case of fast moving particles with arbitrary 
high velocities, such as, for example, particles of ionized gas in some medias at 
very high temperature as : burning reactors, solar winds, nebular galaxies, etc... 
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1) In the non-relativistic case, the original result is due to Carleman, T. , in 
1 ;Di Blasio, G. in 5 Illncr, R. , and Shinbrot, M. in 6 , prove a global result 

in the case of small data and without symmetric assumption ; an analogous 
result is unknown in the full relativistic case. 

2) In the relativistic case, several authors proved local existence theorems, 
considering this equation alone , as Bichteler, K, in [2] ; Bancel, D. ,in [3] , or 
coupling it to other fields equations as Bancel, D. , and Choquet-Bruhat, Y. , 
in g] . Glassey, R. T. , and Strauss, W. , obtained a global result in [7] , in the 
case of data near to that of an equilibrium solution with non-zero density. 

Now it would be interesting to extend to the relativistic case, the global 
existence theorem in the case of small initial data and even for large initial 
data. Mucha, P. B. , studied in [5] and JUj , the relativistic Boltzmann equation 
coupled to Einstein's equations, for a flat Robertson- Walker space-time. 

In this paper, we begin by giving, following Glassey, R. T. , in 8 , the for- 
mulation of the relativistic Boltzmann equation, on a flat Robertson Walker 
space-time background, with a given cosmological expansion factor. We consider 
the homogeneous case which means that the distribution function depends only 
on the time and the momentum of the particles. Such cases are very useful for 
instance in cosmology. Next, we built a suitable functional framework appro- 
priate to the coupling with Einstein's equations, and we construct a sequence of 
operators , to approximate the collision operator, and that provide a sequence of 
real- valued functions that converges to a global solution of the relativistic Boltz- 
mann equation. Moreover, we establish for this sequence, a very simple a priori 
estimation, that leads to a considerable simplification of the method followed 
in [H] . We prove this way, under the only assumption that the cosmological 
expansion factor is bounded away from zero, that the initial valued problem for 
the relativistic Boltzmann equation on the flat Robertson- Walker space-time , 
has an unique global solution, provided that the initial data is sufficiently small, 
and that the solution admits a very simple estimation by the initial data. This 
paper extends the results of [H\ to the case of the curved Robertson- Walker 
background space-time. 

The paper is organized as follows : In section 2, we specify the notations, we 
define the function spaces and we introduce the relativistic Boltzmann equation 
and the collision operator. We end this section by a sketch of the strategy 
adopted. In section 3, we give the properties of the collision operator, some 
preliminary results, and we construct the approximating operators. Section 4 is 
devoted to the global existence theorem. 

2 Notations, function spaces and the relativistic 
Boltzmann equation 

2.1 Notations and function spaces 

A Greek index varies from to 3 and a Latin index from 1 to 3. We consider as 
back-ground space-time, the flat Robertson- Walker space-time denoted (R 4 , g), 
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where, for x = (x a ) = (x°, i') e M 4 , i° = t denotes the time and x = (x l ) the 
space ; g stands for the metric tensor with signature (—,+,+,+) that can be 
written : 

g = -dt 2 + a 2 {t)[{dx 1 f + (dx 2 ) 2 + (dx 3 ) 2 ] (2.1) 

in which, a(t) is a given non-negative regular function of t , called the cosmo- 
logical expansion factor. We suppose a(t) bounded from below, i.e : 

a(t) > Co (2.2) 

where C a > is a given constant. 

We consider the collisional evolution of a given kind of uncharged particles 
in the time oriented curved space-time (R 4 ,<?) . The distribution function of 
the particles we denote by /, is a non- negative real- valued function of both 
the position (x a ) and the 4- momentum p — (p a ) of the particles, and that 
coordinalize the tangent bundle T(R 4 ) i.e : 

/ : T(R 4 ) ~ R 4 x M 4 -> M+, (x a ,p a ) h-> f(x a ,p a ) G M+. 

For p — {p° tP 1 ) — (p° 7 p), and q = (q°,q) we set : 

p.q = a 2 (t)(p 1 q 1 +p 2 q 2 + P 3 q 3 ); \ P | 2 = a 2 {t)((jp x ) 2 + (p 2 ) 2 + (p 3 ) 2 ) (2.3) 

We consider massive particles with a rest mass that can be rescaled torn = 1 . 
The particles are required to move on the future sheet of the mass-shell whose 
equation is g(p,p) = — 1 , or equivalently, using (|2.1Jl and (|2.3(l : 

P° = Vi+TW (2-4) 

Hence, / is in fact defined on the subbundle of T(M 4 ) defined by (2.4) and 
whose local coordinates are x a and p 1 . Now we consider the homogeneous case 
for which / depends only on the time x° = t and p — (p 1 ). The framework we 
will refer to is L 1 (M 3 ) whose norm is denoted || . || ; we set, for r £ R, r > : 

X r = {f€L 1 (R 3 )J>0 a.e,\\f\\<r} 

Endowed with the distance induced by the norm || . ||, X r is a complete and 

connected metric space. Let / be a real interval; we denote by C[I; L 1 (R 3 )] the 

Banach space of continuous and bounded functions g : I — > L 1 (M 3 ) , endowed 

with the norm : ||| g |)|= Sup \\ g(t) \\. We set : 

tei 

C[I; X r ] = {g e C[I; L^M. 3 )} , g(t) e X r , Vt G / } 

Endowed with the distance induced by the norm ||| . ||| ; C[I;X r ] is a complete 
metric space. 
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2.2 The relativistic Boltzmann equation 



The relativistic Boltzmann equation for uncharged particles in the homogeneous 
case we consider , on the Robertson- Walker curved space-time (R 4 ,g) can be 
written : 

with a = 4|, and where Q is the collision operator we now introduce. In the 
case of binary and elastic collisions, Q is defined as follows, p and q standing for 
the momenta of two particles before their collision, p' and q' standing for their 
momenta after the collision, and where / and g are two functions on R 3 : 
1) 

W,.9) = Q + (/,s)-<r (/,<?) (2.6) 

with : 
2) 

Q + (f,9)= [ f ftf)gtf)A<t,p,qj,<ftdu (2.7) 

JR3 1 Js 2 

3) 

Q-(f,9)= f [ f(p)g(q)A(t,p,q,p',q')du J (2.8) 

jr3 q Js 2 

in which; 

4) S 2 is the unit sphere of M 3 whose area element is denoted dcu; 

5) A is a non-negative real valued function of the indicated arguments, called 
the collision kernel or the cross section of the collisions and for which we require 
the boundedness and the symmetry assumptions : 

< A(t, p, q, p', q') < & (2.9) 

A(t,p, q, p', q 1 ) = A(t, p', q', p, q) (2.10) 

where C\ is a given constant. 

6) As consequences of the conservation law p + q = p 1 + q 1 , that gives : 

P° + q =p'° + q' ; p + q=p' + q' (2.11) 
i) we have, using (|2.4|l and the first equation l|2.11|l : 

y/l+ \p\ 2 + y / l+ \q\ 2 = \P'\ 2 + y/l+ \ q' | 2 (2.12) 

which expresses the conservation of the quantity : 

e = \P\ 2 + VU | q \ 2 (2.13) 

called the energy of the unit rest mass particles ; 
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ii) we set, to express the second equation (|2.11|) and following Glassey R.T 

\P' = P+b(p, q,u)uj 

[(? =p-b(p,q,(j)w; ui e S 



(2.14) 



in which b(p, q, w) is a real- valued function. Then, 12.1211 and H2.14fl lead to 
a quadratic equation in b that solves to give the only non trivial solution: 

with p = jjz, e defined by (|2.13() and the dot product by (|2.3() . The Jacobian 
of the change of variables (p, q) i— > (p', (f ) defined by l|2.14|) is given, using the 
usual properties of determinants, by : 

f^) = _^V! (2 . 16) 

d(p,q) p°q° 

Notice that formulae (|2.15|l and l|2.16|l are generalizations to the Robertson- 
Walker space-time, of analogous formulae established in Glassey, R. T. , jS], in 
the case of the Minkowski space-time whose metric corresponds to the particular 
case of <|2.1ll where a(t) = 1. It appears clearly, using (|2~3jl and (j2~T4l) that 
the functions in the integrals (|2.7|l and (|2.8|l depend only on p, q, ui and that 
these integrals give two functions Q + {f,g), Q~(f,g) of the single variable p. In 
practice, we will consider functions / on M 4 that induce for t £ M, functions 
f(t) on M 3 , defined by : /(t)(p) = f(t, p) 

Now, solving the Boltzmann equation (|2.5() that is a first order p.d.e, is 
equivalent to solving the associated characteristic system that writes, taking t 
as parameter: 

f = -2-; f = 4rW,/) (2.i7) 

di a at p° 

We study the initial valued problem for the system (|2.17l) with initial data: 

P l (o) = y l ; /(o) = / (2.18) 

In l|2.17|l . the equations inp 1 integrate directly to give, setting y = (y l ) £ R 3 ; 

P(t,y) = ^-y (2.19) 

whereas the initial valued problem for / is equivalent to the integral equation 
in/: 

f(t,y) = fo(y)+f - Q{Lf){ Sl y)ds (2.20) 

Jo P 

Finally, solving the initial valued problem for the relativistic Boltzmann equa- 
tion 12.5fl with initial data / , is equivalent to solving the integral equation (|2.20() . 
The strategy we adopt to solve l|2.20|l will consist of : 
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a) Constructing an approximating operator Q n , with suitable properties, 
t will converge point-wise in L 1 (R 3 ) to \Q. 

b) Solving by usual methods, the following approximating integral equation: 

f(t,y) = fo(y)+[ Qn(f,f)(s,y)ds (2.21) 



obtained by replacing in (|2.2U|I . -^Q by Q n , and to obtain a global solution /„ 
that will converge to a global solution / of (|2.20|1 . in a suitable function space, 
that will be, given the general framework, a space C[0, +00; X r ] for a convenient 
real number r > 0. Notice that, by (|2.19() . in M. 3 , the volume elements dp, dy 
and the corresponding norms || . \\p and || . \\ y we adopt to be || . ||, are linked 

% 

5 (t) y ' 11 " p a e (t) 



dp=^dy; || . = ^ || . || (2.22) 



3 Preliminary results and approximating oper- 
ators. 

In this section, we give properties of the collision operator Q, we construct and 
we give the properties of a sequence of approximating operators Q n to -^Q . 

Proposition 3.1 If f, g E L\R 3 ) , then Q+(f,g) 1 Q-(f,g) e L^R 3 ) and: 

II^Q + (/,<?)ll<C||/||||<?ll (3.1) 
p 

||^Q-(/,9)||<C||/||||.9|| (3.2) 

II ~ Q + U, f) - - Q + {9,9) ll< CCII / II + II 9 II) II / - 9 II (3-3) 
p° p° 

II ~ Q-{f, f) ~ Q-{9,9) ll< C(|| / || + || g ||) || / - g || (3.4) 

\\- Q(fJ)~- Q(9^)\\<C(\\f\\ + \\ g\\) || f-g\\ (3.5) 
p° p° 

where C = *?£i£±2i 

Proof: 1) the expression l|2.7|l of Q + {f,g) gives, using the upper bound (|2.9|l 
of A : 



-Q + (f,g)\\p<C ia 3 (t) / / \f(p>)\\f(<?)\du; (a) 

P Jr3 y K 3 p g J S 2 
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Now the change of variable (p, q) i — ► (p', q') defined by (|2.14|) and its Jacobian 
(|2.1(j[l give: dpdq — p lo q lo dp'dq' ; then we deduce from (a), using l|2.4l) that gives 

-7^ < 1 : 

p « — 

II IIp< C ia 3 (t) f | /(p') I dp' f | g(q') | dq' [ Aw = 4^a 3 (i) \\ f \\ p \\ 

P° J R 3 J S 2 

(0} then follows from the second relation l|2.22|l and (|2.2|l . □ 

2) The expression ()2.8jl for Q~(f,g) gives, using (|2.9jl and 12.4(1 that gives 

-4» < 1 ; 

II - Q-{L9) \\p< C ia 3 (t) [ | f{p) \dp f | <?(g) | dq [ dw = 47rC ia 3 (0 || / || p -|| 9 

P JR3 Jm 3 JS 2 

and 1)3. 2|l follows from the second relation (|2.22[) and (|2.2|l □ 

3) Inequalities (13.31) . 13. 4|) . I|3.5|l are direct consequences of (|3.1|) . I3.2|l . 

Q = Q + — Q~ , and the fact that, the expressions (|2.7|l and (|2.8|) show that Q + 
and Q~ are bilinear operators, and so we can write, P standing for \Q + or 
jsQ~ ■ P{f, f) - P(9, 9) = P(f, f-g) + P(f ~g,g)- This completes the prove 
of proposition (|3.1|) □ 

We now state and prove the following result on which relies the construction 
of the approximating operators ; C denotes the constant defined in proposi- 
tion (|3~T|) . 

Proposition 3.2 Let r be an arbitrary strictly positive real number. Then, 
there exists an integer n (r) > such that, for every n > n Q (r) and for every 
v £ X r , the equation : 

\friu — — Q(u,u) = v (3-6) 
p° 

has a unique solution u n G X r 

Proof Let r > be given and let v € X r . We can write (3.6), using 
Q = Q+-Q-: 

\/nu — — Q + (u,u) + — Q~(it,u) = v 
p° p° 

Notice that one deduces from the expression (2.8) of Q~ , since p is not a variable 
for that integral, that Q~(f,g) = fQ~(l,g). Hence the above equation can be 
written, using Q~(u,u) = uQ~ (l,u), for u > a.e. 

u = F n (u) (a) 

where 

v+-k;Q+(h,h) 
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(a) shows that a solution u of (3.6) can be considered as a fixed point of the map 
h i — ► F n {h) defined by (b). Let us shows that F n is a contracting map from X r 
to X r . We will then solve the problem by applying the fixed point theorem . 
i)We deduce from (3.1) that F n : X r i — ► L 1 (M 3 ) , and that since v, h € X r : 



\\F n (h) || < 
ii) Let g, h e X r . (b) gives : 

F n (h)-F n (g) = 



^ C || h || || fe || < r(l + CV) 



(c) 



(d) 



1 



>,/>) || 



Since g > 0, ft > a.e., we have , taken n > 0: 

[v^+^Q"(i,^)][v^+^Q"(i,fl)] > «>V^. 

We then deduce from (d), using once more Q~(f,g) = fQ~(l,g), that 

II F n (h) - F n (g) || <|| - o Q+{h,h) - - o Q + {g,g) \\ + \\ - o Q-{v,g) - .. 

+ II ±Q-[±Q + (h,h),g] - ±Q-l±-Q+(g,g), h] || (e) 

We can write, using the bilincarity of Q~ : 

~o Q~( v > 9) -j Q-{v,h) = j o Q-{v,g-h) 

So, in (e) we apply (3.3) to the first term, (3.2) to the second term to obtain, 
using || g \\< r, || h \\< r, \\ v \\< r : 

|| ±.Q+(h,h)~Q+(g,g) ll + H - o Q-{v,g) + - o Q-{v,h) \\<3Cr\\h-g (f) 

pO pO pU pU 

Concerning the last term in (e), we can write, using the bilinearity of Q~ : 



Q- 



Q+(h,h),g 



pO 



■Q + (g,g),h 



l 



+ -Q- 

po 



Q+(h 7 h),g-h 



l-Q+(h,h)-^Q+(g,g),h 



We then apply to the right hand side of the above equality : 

1) For the first term, property (3.2) of Q~ , followed by property (3.1) of Q + 

2) For the second term, property (3.2) of Q~ , followed by property (3.3) of 
Q + to obtain, using || g \\< r, \\ h \\< r : 



-Q + (h,h),g 



<SC^\\h-g\\ (h) 
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thus (e) gives, using (f) and (h) : 

II F n (h)-F n (g) ||< 3Cr(1 J i Cr) \\h-g\\ (i) 

It is then easy , using (c) and (i), to choose an interger n (r) such that, for 
n > n (r) , we have: 

1 + Cr , 3CY(1 + Cr) 1 

=- < 1 and v — '- < - (3.7) 

y/n ~ " 2 

This shows , using once more (c) and (i) that, for n > n (r), F n is a contracting 
map of the complete metric space X r into itself. By the fixed point theorem, 
F n as a unique fixed point u n e X r that is a unique solution of (3.6). 
In all the following n (r) stands for the integer defined above. 
The above result leads us to the following definition : 

Definition 3.1 Let n G N, n > n (r) be given. 

1) Define the operator R(n,Q) : X r — > X r , u i — > R(n,Q)u as follows : 
for u G X r , R(n,Q)u is the unique element of X r such that : 

VnR(n, Q)u ~ —Q [R(n, Q)u, R(n, Q)u] = u (3.8) 

pO 

2) Define the operator Q n on X r by : 

Q n (u,u) — ny/nR(n,Q)u — nu (3.9) 
We give the properties of the operators R(n, Q) and Q n . 

Proposition 3.3 : Let r £ R^ 7 n, m G N m > n Q (r), n > n D (r), be given. Then 
we have, for every u, v e X r 

|| Vn~R(n,Q)u || = || u \\ (3.10) 

II y/nR(n, Q)u - u II < — (3.11) 

n 

|| R{n,Q)u - R(n,Q)v \\ < || u - v \\ (3.12) 

V n 

Q n {u,u) = ^Q[Vn~R(n,Q)u,VnR(n,Q)u] (3.13) 

|| Q n {u,u) - Q n {v,v) \\ < K \\ u - v || (3.14) 

|| Q n (u,u)- \Q(v,v) I) < -+K\\u-v\\ (3.15) 
p° n 

K K 

\\ Q n (u,u) - Q m (v,v) \\ < K\\u-v\\^ 1 (3.16) 

n m 

Here K = K{r) is a continuous function of r. 
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Proof. 1) (3.10) will be a consequence of the relation: 

/ - Q{f,9)dp = Q, V/^gei^M 3 ) (3.17) 

we now establish. This is where we need assumption (2.10) on the kernel A. We 
have, using (??), definition (2.7) of Q + , the change of variables (p, q) i — > (p'q') 
defined by (2.14) and its Jacobian (2.16), and since p° = y/T+ \ p | 2 : 

W(f,9)dp = a*(t) [ dp' [ dq' [ ff)9ti)S&,?,p,q) 

3 P° V R 3 P U J S 2 ^l+|p'|Vl+U'| 2 



On the other hand, we have, using definition (2.8) of Q and p° = yl+ | p | 2 : 

W ±- Q-(f, 9 )dp=a\t)[ dpi dof f j^^f;%^ 
p° 7m3 Jk3 7 S 2 v/i+ I p | 2 ^i+ I Q | 2 

It then appears that I = J and 13.17|l follows from Q = Q + — Q~ . Then, 
to obtain (3.10), we integrate (3.8) over M 3 with respect to p, use (3.17), and 
conclude by (2.22) that gives the equivalence between || . \\ p and | . j| □ 

2) The definition (3.8) of R(n, Q) gives, using property (3.5) of Q with 
/ = R(n, Q)u and g = : 

|| ^R(n,Q)u-u ||=|| —Q[R{n,Q)u,R(n,Q)u] \\<C\\ R{n,Q)u \\ 2 
p° 

and (3.11) follows from (3.10) that gives || R{n,Q)u ||< M < 

3) Subtracting the relations (3.8) written for u and v gives: 

u — v \ i 
i?(n, Q)u—R(n, Q)v — — -j= — V—j= — Q [R(n, Q)u, R(n, Q)u] y= — Q [R{n, Q)v, R(n, Q)v] 

which gives, using property (3.5) of Q and once more (3.10): 

II ii — v II ICr 

|| R(n, Q)u - R(n, Q)v \\< ^ + -=- || R(n, Q)u - R(n, Q)v \\ 

\Jn y/n 

(3.12) then follows from n > n {r) and (3.7) that implies < \. 

4) The expression (3.13) of Q n is obtained by multiplying (3.8) by n, using its 
definition (3.9) and the bilinearity of Q. 

5) Property (3.14) is obtained by using (3.13), property (3.5) of Q, followed by 
properties (3.10) and (3.12) of R(n,Q). 

6) Notice that by (3.10), y^riR(n, Q)u 6 X r since u £ X r . Then, the expression 

(3.13) of Q n and property (3.5) of Q give : 

|| Q n (u,u)-^-Q{v,v) \\< 2Cr \\ ^R{n,Q)u-v \\< 2Cr (|| sfriR{n, Q)u - u \\ + \\ u-v \\) 
p° 
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hence, (3.15) follows from (3.11). 
7) The expression (3.13) of Q n gives: 

Q n {u, u)—Q m (v, v) = -^Q [y/nR(n, Q)u, y/nR(n, Q)u] — -^Q [y/mR(m, Q)v, y/mR(m, Q)v] 

then, (3.16) follows from (3.5), (3.10) and (3.11), adding and subtracting u and 
v. □ 

Remark 3.1 . (3.15) shows by taking u = v, that Q n converges pointwise in 
the L 1 (M. 3 )-norm to the operator ^Q- From there, the qualification of 
"approximating operators" we give to the sequence Q n . 

We now have all the tools we need to prove the global existence theorem. 



4 The global existence theorem 

With the approximating operator Q n defined by (3.8), we will first state and 
prove a global existence theorem for the "approximating" integral equation: 

l-ta+t 

f(t + t,y)=g (y)+ Qn(f,f)(s,y)dy, t > 0, (4.1) 



in which t is an arbitrary positive real number, and g a is a given function that 
stands for the initial data at time t Q \ for practical reasons,we study l|4.1|l rather 
than (|2.21(l that is the particular case when t a = and g a — f a . Next, we prove, 
by a convenient choice of t that the global solution /„ of 1)4.1(1 converges in a 
sense to specify, to a global solution / of 1(2.20(1 . 

Proposition 4.1 Let r be an arbitrary strictly positive real number. Let g a € 
X r and let n be an integer such that n > n a (r). Then, for every t a €E [0, +oo[, 
the integral equation j4-l\ ) has a unique solution f n € C[t a , +oo; X r ]. Moreover, 
f n satisfies the inequality : 

III /n|||<|| 3o|| (4.2) 

Proof: 

We give the proof in 2 steps. 

Step l:Local existence and estimation. 

Consider equation ((4.1(1 for t 6 [t ,t + 5] where 5 > is given. One verifies 
easily, using (3.8) that gives : Q n {f, /) + nf = n 2 R(n, Q)f and 

f t [e nt f(t + t, y)] =e nt f t + nf (t a + t, y), that (gjj is equivalent to : 

f(t + t,y) = e- nt g (y)+ [ W^-"^ R(n,Q)f(t a + s,y) dy (4.3) 



We solve 1(4.3(1 by the fixed point theorem. Let us define the operator A over 
the complete metric space C[t , t Q + 5; X r ) by : 



Af(t + t,y) = e- nt 9o (y) + / n^e-"^ R(n, Q)f(t + s, y) dy (4.4) 
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i) Let / e C[t ,t + S;X r ]; 14.4|1 gives, since || g \\< r and using (3.9) that 
gives, || ^/nR(n, Q)f(t + s) \\ = \\ f(t a + s) \\ < r, and for every t £ [0, 8} : 

II Af(t a + t)\\< e- nt r + nre~ nt f e ns ds = e~ nt r + r e - nt {e nt - 1) = r 

Jo 

Hence 1 1| Af \\\< r, and this shows that A maps C[t a , t a + 8; X r ] into itself. 

ii) Let f,geC[t o ,t o + S;X r ];t€[0,6}; (4.4)gives: 

(Af - Ag) (t +t, y) = [ nV^e~ n ^ [R(n, Q)f - R(n, Q)g] (t a +s, y) dy 
Jo 

which gives, using property (3.12) of R(n,Q) and e -n '* -s ) < 1 : 
\\\Af-Ag\\\<2n5\\\f-g\\\. 

It then appears that A is a contracting map in a metric space C[t ,t + S;X r ] 
where 2nS < i, i.e. 8 S]0, JM. Taking 8 = 4-, we conclude that A has a 
unique fixed point f° £ C[t ,t + j-;X r ], that is the unique solution of (4.3) 
and hence, of (4.1). 

Now we have, by taking t — in (4.3), f(t ,y) = 9o{y)- The solution f° then 
satisfies: 

f°{t +t,y) = e- nt f°(t , y) + [ n^e~ n{t -^ R(n, Q)f°{t + s, y) dy 

Jo 

which gives, multiplying by e nt , using once more (3.10), and for t G]0, y-]: 

II e nt f°{t + t)\\< || /» || +n f || e ns fZ(t + s) \\ ds 

Jo 

this gives, by Gronwall Lemma: e nt \\ f°(t +t) \\< e nt \\ f°(t ) |j and hence 

III fn III < II fnito) || (4.5) 

Step 2: Global existence; estimation and uniqueness. 
Let k £ N. By replacing in the integral equation (4.1) t by : 

*° + in ' *° + in ' ' ' ' ' *° in ' " ' ' ste P ^ te ^ s us *^at 011 eac ' 1 m t erva l 

If. = [t + j-, t Q + j- + j-] whose length is the initial value problem for the 

corresponding integral equation has a unique solution /* G C[Ik,X r ], provided 

that the initial data we denote /* \t a + J-] is a given element of X r ; /* then 

satisfies: 

/n(*o + h + *> ?/) = /n(*o + y) + Jj+j: + * Qn (/£, /*) («, 2/) ^ 

and (4. 5) implies: 

|||/,tlll<ll^(< + ^) II- 
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Now, [0, +oo[= U [to + £ , t + £ + t] ; Define : 

feeN 

/„ fc (io + ^,y) = /r i (to + ^,y) if fc>land f° n {t )=g . 

Then, the solution J* -1 and /* that are defined on [t + , i Q + ^] and 
[to + *o + ^] overlap at t = t„ + 
Define the function /„ on [0, +00 [ by 



/«(*) = /£(*) if te 



fc fc + 1 



Then, a straightforward calculation shows, using the above relations for 

k, k — 1, k — 2, • • • ,1,0, that /„ is a global solution of (4.1) on [0, +oo[ satisfying 

the estimation (4.2) and hence, /„ e C[t , +00; X r \. 

Now suppose that /, g € C[t , +00; X r ] arc 2 solutions of (4.1), then: 



if -9) (to + t,y) 



I [Qn(f,f)-Qn(g,g)}(t + s,y)ds i>0. 

This gives, using property (3.14) of Q n : 

\\{f -g){to + t)\\< K ( \\(f-g)( to + s)\\ds t>0 

Jo 

which implies, using Gronwall Lemma, f = g, and the uniqueness follows. This 
ends the proof of proposition 4.1.D 

We now state and prove the global existence theorem. 

Theorem 4.1 Let f Q G L 1 (R 3 ), f a > 0, a.e; be arbitrary given. Then the 
Cauchy problem for the relativistic Boltzmann equation on the Robertson- Walker 
space-time, with initial data f , has a unique global solution f € C[0, +00; i 1 (M 3 )] 
f{t) > 0, a.e,Vt G [0, +oo[; / satisfies the estimation: 

Sup || /(t) || < || /„|| (4.6) 
te[o,+oo[ 

Proof: Fix a real number r >|| f \\. We give the proof in 2 steps. 
Step 1: Local existence and estimation. 

Let t > and g e X r be given. The proposition 4.1 gives for every integer 
n > n (r), the existence of a solution /„ e C[t , +00; X r ] of (4.1) that satisfies 
the estimation (4.2). It is important to notice that this solution depends on n. 
Let T > be given; we have /„ € C[t ,t + T;X r ], Vn E N, n > n (r). Let 
us prove that the sequence /„ converges in C[t , t + T; X r ] to a solution of the 
integral equation : 

/(t„ + t, y) = 9o {y) + f ° + - Q Q (/, /) (s, y) dy, t e [0, T] : (4.7) 
Jt„ P 
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Consider 2 integers m,n > n (r) ; we deduce from (4.1), that, Vt € [0, T]; 

/■t +< 

/n(io + t,y)- f m (t +t,y)= / [Q n (/„, /„) - Qm (/m, /m)] («, y) dy 



this gives, using the property (3.16) of the operators Q n , Q m and Gronwall 
lemma: 1 1| /„ — f m ||| < + ^) KTe KT . This proves that /„ is a Cauchy 
sequence in the complete metric space 

C[t a , t Q + T, X r ). Then, there exists / g C[t a , t Q + T, X r ] such that: 

/„ converges to / in C[t a , t a + T, X r ] (4.8) 
The definition of X r implies that f(t a + t) > a.e; 

Let us show that / satisfies (4.7). The convergence (4.8) implies that, f n (t + t) 
converges to f(t +t), Vt g [0, T]; then f n (t ) converges to f(t ); but f n (t ) = g a , 
then f(t ) = g ; next we have, using property (3.14) of Q n and Q : 



r-t +t 

J to 



Qn(fnJa)--Q(.f,.f) 
Po 



(s,y)dy\\< ^1+KT\\\ /„-/ 



it. 

which shows that 

rto+t rto+t 

J Q„ (/„,/„ )(s)ds converges to J Q(f,f)(s)ds in L 1 (R 3 ). 

Hence, since /„ satisfies (4.1), / satisfies (4.7), Vt <E [0,T] and a.e. in y g M 3 . 
Finally, (4.2) and (4.8) imply: 

III / III < II /(*o) II (4-9) 
Step 2: Global existence, estimation and uniqueness 

Since t € [0, +oo[ is arbitrary and since Vn > n (r), the solution f n of (4.1) 
is global on [t a , +oo[, Step 1 tells us that, given T > 0, by taking in the integral 
equation (4.7), t Q = 0, t a = T, t a = 2T, • • • , t Q = kT, ■ ■ ■ k g N, then, on each 
interval Jfe = [kT, (k + 1)T], fc g N, whose length is T, the initial value problem 
for the corresponding integral equation has a solution f k e C[Jk, X r ], provided 
that the initial data we denote f k (kT), is a given element of X r . 
We then proceed exactly as in Step 2 of the proof of Proposition 4.1 by writing 

for T > given, that [0, +oo[= (J [kT, (k + 1)T], to overlap the local solutions 

feeN 

f k g C[kT, [k+l)T;X r ], by setting f k (kT) = f k ~\kT) if k > 1 and /°(0) = / , 
and we obtain a global solution / e C[0, +oo; X r ] of the integral equation 
(2.20) defined by f(t) = f k (t), if i g [kT,(k + 1)T]. Next, (4.9) that gives: 
III / fc III < II f k (kT) ||, shows that / satisfies the estimation (4.6). 
Now if /, g e C[0, +oo; i 1 (R 3 )] are 2 solutions of (2.20) for some initial data f a , 
(4.6) and property (3.5) of Q give: 

||(/-s)(t)||<2C||/ || [ t \\(f-g)(a)\\ ds, Vt>0. 

Jo 
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which gives, by Gronwall Lemma, / — g, and the uniqueness follows. 
Now as we indicated, solving the initial value problem for the relativistic Boltz- 
mann equation is equivalent to solving the integral equation (2.20). This ends 
the proof of theorem 4. ID. 

Remark 4.1:In a future paper, we study the Boltzmann equation coupled to 
the Einstein equations in the case of several species of particles with different 
rest masses. 
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